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找 性 航天 器 非 约束 模 态 动力 学 模型 的 降 阶 研究 


(北京 信息 科技 大 学 机 电工 程 学 院 ,100192 北京 ) 


摘 要 :针对 中 心 刚 体 - 挠 性 附件 的 挠 性 航天 器 ,应 用 拉 格 朗 日 方程 推导 出 动力 学 模型 ,再 利用 模 态 
展开 法 得 到 非 约束 模 态 动 力学 方程 ,然后 分 别 基 于 模 态 价值 分 析 准 则 和 内 平衡 降 阶 准则 的 相关 理 
论 探讨 模型 降 阶 。 通 过 模 态 价值 分 析 准 则 能 够 显示 出 各 阶 模 态 占 系统 总 价值 的 份额 ,保留 贡献 较 
多 的 模 态 组 成 降 阶 系统 。 内 平衡 降 阶 准则 通过 非 奇 异 线性 变换 ,将 系统 模型 转换 为 内 平衡 系统 ,在 
内 平衡 系统 内 会 掉 Gram 短 阵 对 角 线 上 较 小 的 值 所 对 应 的 模 态 以 构成 降 阶 系统 。 最 后 ,利用 数值 仿 
真 对 中 心 刚体 - 单 侧 挠 性 梁 航 天 器 的 非 约束 模 态 动 力学 模型 进行 模型 降 阶 ,通过 两 种 准则 对 系统 模 
型 降 阶 数值 计算 ,结果 表明 两 种 降 阶 准则 均 能 够 有 效 地 对 原 系 统 进行 模型 降 阶 。 
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Modal reduction analysis of unconstrained modal dynamic 
of large flexible spacecraft 


WEI Zhidong, GE Xinsheng 


(Mechanical Electrical Engineering School, Beijing Information Science & Technology University ,100192 Beijing,China) 


DAbstract :A dynamical model was derived by using Lagrange equation for a spacecraft with a center rigid 
body and a single large flexible structure ,and then the unconstrained modal dynamic equation was devel- 
:oped by modal expansion. The relevant theoretical knowledge of modal cost analysis criterion and internal 
balance technique was introduced to solve model reduction. Modal cost analysis approximation formula to 
solve a weak damping and sufficiently separated from natural frequency in system was given in this paper. 
Modal cost analysis criterion could show the share of each modal in the total value of the system ,and the 
reduced order system was composed of the modes with more contributions. Internal balance technique 
meaning the state space system was converted to internal balance system by non-singular linear transforma- 
tion. The reduced order system was constructed according to discard the modes corresponding to the values 
on the diagonal Cram matrixes of the internal balance system. Finally ,numerical simulation was applied to 
reduce the order of the unconstrained modal dynamics model of the spacecraft with central rigid body and 
a single flexible beam. The numerical calculation of the reduction of the system model by two criteria shows 


that the original system is effectively reduced. 
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随 着 航天 科技 不 断 发 展 ,大 挠 性 结构 航天 器 动 
力学 模型 "31 成 为 大 家 普遍 关注 的 问题 。 早 期 挠 性 
航天 器 大 多 都 是 中 心 刚体 的 惯量 比 挠 性 附件 惯量 大 
许多 , 建 模 时 通常 将 中 心 刚体 视 为 静止 不 动 。 当 搁 
性 航天 器 中 心 刚体 的 惯量 与 挠 性 结构 的 惯量 比值 相 
接近 时 , 挠 性 结构 如 果 继续 采用 悬臂 状态 的 传统 方 
法 对 大 挠 性 结构 的 航天 器 进行 动力 学 分 析 , 将 会 产 
生 偏差 中。 针对 大 挠 性 航天 器 , Wang 等 中 采用 有 
限 元 法 及 拉 格 朗 日 法 建立 了 大 挠 性 航天 器 的 姿态 动 
力学 方程 。 周 志 成 等 采用 有 限 元 法 描述 了 柔性 辟 
杆 和 太阳 沟 , 通 过 拉 格 朗 日 法 推导 了 柔性 组 合 航天 
器 本 动力 学 方程 。 非 约束 模 态 定义 为 中 心 刚体 不 固 
矫 整 个 航天 器 无 阻尼 的 自由 振动 过 程 。 文 献 [6-8] 
据 时 了 非 约 束 模 态 动量 系数 的 恒等式 关系 ,并 给 出 
拨 性 航天 器 非 约束 模 态 动力 学 方程 。Hablani'” | 描 
述 巴 非 约束 模 态 为 不 存在 外 力 和 阻尼 的 情况 下 , 航 
苍 器 的 自由 振动 过 程 , 且 推 导 了 非 约 束 模 态 的 正 交 
性 昌 旺 等 "给 出 了 拨 性 航天 器 在 轨 的 非 约束 动 
志江 模型 ,得 到 了 非 约束 模 态 频率 般 大 于 约束 模 
埠 颈 率 的 结果 ,其 原因 是 与 约束 模 态 相 比 , 非 约束 模 
态 栅 当 于 降低 了 约束 ,使 找 性 附件 振动 更 加 自由 , 因 
此 颖 率 升 高 ,特别 是 在 模 态 基 频 上 比较 明显 。 文 献 
[E12] 采 用 哈密 顿 原理 推导 了 拨 性 航天 器 动力 学 
方程 ,并 给 出 了 航天 器 系统 非 约束 模 态 振 型 的 正 交 
性 和 求解 非 约束 模 态 振 型 的 方法 。 用 约束 模 态 法 和 
非 的 束 模 态 法 中 对 比 研究 了 搁 性 航天 器 建 模 相关 
问题 ,并 讨论 了 约束 模 态 频 率 和 非 约 束 模 态 频率 的 
计算 误差 相关 问题 。 

针对 挠 性 结构 航天 器 ,由 于 挠 性 结构 的 引入 使 
得 动力 学 模型 的 自由 度 变 得 很 大 ,而 主动 控制 设计 
的 实现 则 需要 系统 模型 的 阶 数 尽 可 能 的 低 ,这 就 需 
要 对 模型 进行 降 阶 处 理 ,用 低 阶 模型 近似 代替 原 高 
阶 模型 。 模 型 降 阶 方法 通常 采用 以 下 3 种 降 阶 准 
则 :惯性 完备 性 准则 "中 根据 约束 模 态 惯性 完备 性 
恒等式 计算 各 阶 模 态 的 模 态 质量 和 转动 惯量 ,通过 
忽略 各 阶 模 态 质量 和 惯量 对 系统 影响 较 小 的 模 态 ， 
从 而 达到 模 态 截断 效果 ,该 方法 理论 发 展 成 熟 、 降 阶 
较为 方便 ,常用 于 建 模 初 期 ; 模 态 价值 分 析 准 则 "1 
需要 计算 系统 各 阶 模 态 的 模 态 价值 ,各 阶 模 态 的 模 
态 价值 大 小 一 般 有 数量 级 的 差异 ,只 有 其 中 一 两 个 


模 态 起 重要 作用 ,该 方法 物理 意义 明确 ,使 用 范围 
广 ,形式 简洁 ,是 与 控制 目标 结合 起 来 一 种 广泛 用 于 
各 种 挠 性 结构 的 降 阶 方法 ,其 降 阶 效 果 已 得 到 了 验 
证 ; 内 平衡 降 阶 准则 王 2 是 通过 将 动力 学 模型 转换 
到 状态 空间 系统 ,再 将 状态 空间 系统 转换 为 内 平衡 
系统 ,在 该 系统 内 通过 舍 掉 Cram 对 角 阵 对 角 线 上 较 
小 的 值 所 对 应 的 模 态 ,从 而 达到 模 态 截断 效果 。 

本 研究 针对 非 约 束 模 态 动力 学 模型 进行 降价 分 
析 , 计 算 其 非 约 束 模 态 频率 ,应 用 模 态 价值 降 阶 准则 
和 内 平衡 降 阶 准则 对 非 约束 模 态 动力 学 模型 进行 模 
型 降 阶 。 通 过 对 非 约 束 模 态 动力 学 模型 降 阶 数 值 仿 
真 ,仿真 实验 表明 无 论 是 模 态 价值 降 阶 准则 或 内 平 
衡 降 阶 准则 与 原 系统 相 比 较 都 能 起 到 了 较 好 的 降 阶 
效果 。 


1 动力 学 方程 


考察 中 心 刚体 带 单 侧 搁 性 梁 ,如 图 1 所 示 , 建 立 
惯性 系 oxyz 和 附件 坐标 系 0;x;y;z;,0z 和 oiz; 垂直 纸 
面向 外 。oixiyi2; 的 原点 在 中 心 刚 体 和 挠 性 深 的 连接 
点 o 上 ,oxyz 的 原点 在 系统 的 质心 。 上 ,假设 oy 轴 
方 加 的 位 移 为 y, , 绕 0z 轴 的 转动 角度 为 0, ,附件 坐 
标 系 下 挠 性 梁 上 任意 一 点 4 在 ox; 轴 上 的 坐标 为 
i, 梁 的 搁 度 为 w(xi,t) ,点 4 在 ox 轴 上 的 坐标 为 
x0o。 系 统 质心 o 到 中 心 刚体 和 挠 性 梁 的 连接 点 o; 的 
距离 为 m ,质心 C 为 中 心 刚体 的 质心 。 本 研究 仅 考 
虚 搁 性 梁 的 横 问 振动 ,不 考虑 其 扭转 振动 与 伸缩 
振动 。 

在 惯性 系 下 梁 的 挠 度 w(x; ,六 为 


0,xoro 
Ww(Xi,t) =y, +xo0, + 1 
人 ) 4 " a ,Xo > To ( 


挠 性 梁 任 意 一 点 的 搁 度 用 模 态 展开 为 
DAE )m(t) = Bx) nt) (2) 


其 中 :g(x ) 为 摘 性 梁 的 第 上 阶 非 约 束 模 态 阵型 ; 
人 (站 为 搁 性 梁 的 第 阶 模 态 坐 标 。 
在 惯性 系 下 的 速度 表示 为 


W(x; ,t) = 


0,xoro 


0 = =7raobt| (3) 


w(xi,t) ,Xo >70 
系统 动能 可 写 为 
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T=7T,+T, (4) 了 = 了 -了 为 Lagrange 郧 数 ; 设 广义 力 R 的 分 力 FT、 
其 中 了 表示 为 中 心 刚体 的 动能 。 f 分 别 为 外 力 、 外 力矩 和 模 态 力 ,R=[F 7 f]"。 
也 = (7 +xo0, ) (7， + xoO, ) dm, (5) 将 式 (7) OT 4 
my,+c 0,.+pn =F 
其 中 了 表示 为 挠 性 梁 的 动能 _.. m 
cya+J0,. +hr = 了 (10) 
7.= 记 | (7 +x00, tw (x,t) ) 3 Ap tf 
(ya txo0, + w(xi,t) ) dm (6) 将 方程 式 (10) 写成 矩阵 形式 , 即 
Mai +Kg=R (11) 


图 1 中 心 刚体 - 单 侧 找 性 梁 结 构 


Fig.1 The structure of center rigid body 


and a side flexible beam 


CN 将 式 (5) ~ (6) 代 入 式 (4) ,可 得 
QI7T= 


Fm Yo + Os+yp N+Ohn t+ 

1 二 

FI N+E YO (7) 
HB:m 为 系统 总 质量 , 且 m =m, +mii7 为 找 性 梁 


的 铝 态 坐标 列 阵 ;p 为 非 约束 模 态 动量 系数 , 且 
p 志 | Bdmi;h 为 韭 约束 模 态 角 动 量 系 数 , 且 


hh= | #6q@Bdmi;J 为 相对 于 系统 质心 。 的 转动 惯量 
表示 为 J=J +J;5 5 分 别 为 挠 性 梁 和 中 心 刚体 相 
对 质心 C 的 静 矩 , 且 5 = 2 + +5 = | iodm + | 5odm:= 


| tdmo 
航天 器 势 能 为 找 性 梁 的 变形 势能 , 即 


y= B(x) KB(x) = (8) 
其 中 4 为 挠 性 染 固 有 频率 组 成 的 对 角 阵 。 
根据 Lagrange 方程 
fol) aL) 
Us 9) 
其 中 :4 为 广义 坐标 列 向 量 , 且 g=[y，0，。，1]'; 
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2 模型 降 阶 准则 


模 态 价值 分 析 准 则 和 内 平衡 降 阶 准则 在 各 种 航 
天 器 挠 性 结构 模型 降 阶 研究 中 已 得 到 了 验证 。 现 简 
要 介绍 两 种 准则 的 相关 理论 。 

将 动力 学 方程 式 (11) 的 响应 输出 方程 以 位 移 
坐标 和 速率 表示 为 


y=Ciq+C.4 
其 中 C,\C, 分 别 为 位 移 ,速率 观测 矩阵 。 


对 方程 式 (11) 添 加 阻尼 项 后 与 方程 式 (12) 联 
可 得 系统 广义 的 动力 学 方程 , 即 


0 


(12) 


. (13) 

y=Ciq+C.g 
其 中 MD、K 为 系统 质量 阻尼、 刚度 和 矩阵, 且 M = 
M'" >0,D=D'=0,K=K 0;u 为 系统 的 输入 或 外 


界 干扰 ;B 为 u 的 分 布 矩 阵 。 
2.1 模 态 价值 分 析 准 则 


模 态 价值 分 析 是 通过 计算 系统 中 各 阶 模 态 的 模 
态 价值 和 各 阶 模 态 价值 占 系 统 总 价值 的 比例 的 多 少 
来 对 模 态 的 取舍 ,保留 贡献 较 多 的 模 态 组 成 降 阶 系 
统 ,从 而 实现 对 系统 进行 模型 降 阶 。 

采用 正则 模 态 矩阵 了 ,将 方程 式 (13 ) 转换 到 模 
态 坐 标 空间 , 即 


bee 


(14) 
》 于 Can + Cn 


其 中 :wo 纪 分 别 为 第 i 阶 非 约束 模 态 频率 和 阻尼 


] 力 学 学 报 


第 6 其 


比 ;B=T'B.C, =CJT.C, = C.T。 青 将 方程 式 (14) 
转换 到 状态 空间 , 即 


， =Ax +Bu (15) 
y=Cx 
其 中 , = (ni ”CO1711 0 
-206;: -ow; T 
4=diag(4),4=| 0 |:8=[B!, 


…,B'],B,=[b，0] ,5b, 是 B 的 第 i 行 ;C=[C， 
Cau/ wi y ;CG , Ca/ ,| [e) 

由 方程 式 (15) 可 知 ,系统 的 第 i 阶 模 态 价值 数 
值 解 为 


V,=tr [XC QC]; (16) 
其 中 8 为 权 和 矩阵 ;为 状态 协 方差 矩阵 ,其 满足 Lya- 
pinov 方程 的 解 , 即 
C XA' +AX+BB' =0 (17) 
CG) 系 统 总 价值 了 为 各 阶 模 态 价值 之 和 , 即 
J V = > V, (18) 
CCV 各 阶 模 态 价值 占 系统 总 价值 的 比例 为 
| (MOT) ;= V/V (19) 


2 人 4 内 平衡 降 阶 准则 
ON 内 平衡 准则 实际 上 是 将 状态 空间 系统 变换 到 内 
平衡 系统 ,使 得 内 平衡 系统 的 可 控 可 观 Gram 矩阵 具 
部 对 角 阵 结构 ,从 而 保留 该 系统 内 可 控 可 观 Gram 矩 
阵 对 角 线 上 较 大 的 值 ( 即 奇异 值 o,) 对 应 的 模 态 以 
构 咸 降 阶 系统 ,从 而 实现 对 系统 进行 模型 降 阶 。 
于 方程 式 (15) 的 可 控 和 可 观 Gram 矩阵 分 别 为 
W.W, ,可 由 Lyapunov 方程 得 2 
AW.+WA' +BB' =0 
pe +WA+CC'=0 
设 非 奇异 线性 变换 , 即 
VT (21) 
将 方程 式 (21) 代 入 方程 式 (15 ) 系统 ,可 得 内 平 
衡 系 统 为 


(20) 


ti-Atz+Bu 
(22) 
y=Cx 
其 中 4 =7T7A47T,,B=7T-'B,C =CT,，。 
此 时 该 系统 Gram 矩阵 本 、W, 满足 
WwW. = WwW. =diagle;} (23 ) 
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由 方程 式 (20) 可 知 ,可 控 Gram 矩阵 W. 满 
足 52] 
4 本 + WA +BB'=0 (24) 
将 式 (24 ) 左 乘 和 矩阵 了, 右 乘 矩阵 T… ,同时 将 
A =7T-'AT,,B =7T-1B,C =CT, 代入 ,可 得 
A(T,W.T') + (TW.T')A +BB'=0 (25) 


令 W. = TW.TI, 即 W. = Ti'W. (Ti')", 同 理 


可 导出 可 观 Gram 矩阵 , 即 
W, =TiW.T, (26) 

将 可 观 可 控 Gram 和 矩阵 WW. 、W, 相 乘 ,得 
W.W, =Tr" (W.W,)T, (27) 


可 知 ,和 矩阵 T, 是 矩阵 WW.W, 的 特征 向 量 , 系统 
的 特征 值 为 Hankel 奇异 值 0 。 


3 数值 仿真 


如 图 1 所 示 ,考察 中 心 刚体 带 单 侧 挠 性 粱 ,中 心 
刚体 为 盒 型 立方 体 , 边 长 为 1.0 m, 其 他 材料 参数 
和 挠 性 粱 上 … 的 材料 参数 如 表 1 所 示 。 

表 1 中 心 刚 体 和 挠 性 梁 材 料 参数 
Tab.1 Material parameters of center rigid 


body and flexible beam 


中 心 刚 体 挠 性 梁 
参数 值 参数 值 
密度 /(kg* m) 7 800 密度 /(kg .mm 一) 920 
质量 /kg 231.67 质量 /kg 5.21 
弹性 模 量 /Pa 2.0 x10!! 弹性 模 量 /Pa 2.6x10! 
泊 松 比 0.3 泊 松 比 0.3 
惯量 / diag 惯量 / diag 
(kg .m2) 163.71 63.71 63.71| (kg .m2) |40.940.42 41.36| 
刚体 壁 厚 /m 0.005 长 / 宽 / 厚 /m 4/1/0.015 


分 别 采用 模 态 价 值 分 析 准 则 和 内 平衡 准则 验证 
挠 性 航天 器 模型 降 阶 的 有 效 性 和 可 行 性 。 假 设 挠 性 
梁 的 自由 端 端点 分 别 受到 脉冲 激励 和 正弦 激励 响 
应 。 梁 端点 的 初始 位 移 为 0.01 m, 初 始 速 度 为 0, 撤 
除 激励 后 挠 性 梁 做 自由 衰减 振动 。 分 别 求 得 非 约束 
模 态 平 动 耦合 系数 和 转动 耦合 系数 ( 表 2 ~3) 、 找 性 
梁 的 约束 模 态 频率 w, 和 由 方程 式 (11) 计 算 的 系统 
非 约束 模 态 频率 wx( 表 4) 。 
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表 2 系统 前 10 阶 非 约 束 模 态 平 动 看 合 系数 p, 


Tab.2 Translational coupling coefficient p;, of unconstrained 


mode of the lst to 10th of flexible spacecraft system 


j plj,x) /kg p(j,y) 《kg p(j,z) /kg 
1 -0.5998x107? 0.1831x10™3 1.1300 
2 -0.2034x107™3 0.2864x10-? -0.682 1 
3 0.2227 x10-3 -0.8948x10-4 —0.0338 
4 -0.6334x10-5 -0.1202x10-? 0.5148 

5 -0.3888x10-3 -0.1177x10-4 -0.701 1 x107? 
6 0.2719 x10-4 0.7202 x10-3 -0.2006 
7 0.1806x10-; 0.2071x10- 0.013 8 

8 -0.4859x10-4 -0.4874x10-3 0.285 9 
9 -0.1121x10-3 -0.1845x10- 0.019 4 
IO 0.9671x10-4 0.342 8 x10-3 -0.1100 


J 表 3 ”系统 前 10 阶 非 约束 模 态 转动 厅 合 系数 有 


CY. 3 Rotational coupling coefficient h;, of unc i 


mode of the lst to 10th of flexible spacecraft system 


h(j,x)/ h(j,y)/ h(j,z)/ 
(kg m )'® (kg om (kg om) 
3.7440 0.0109 0.2535 x107? 
-1.1870 -0.042 8 0.3660x10-? 
0.080 1 0.332 2 0.6601x107? 
1.0660 0.0160 -0.1639x107? 
0.7253 x107? -0.2145 -0.5316x10-3 
-0.4050 -0.0190 0.593 1 x10-3 
-0.1028 0.0308 -0.7057x10-4 
0. 8067 -0.8881x10-” -0.8643x10-3 
0. 1687 -0.0605 -0.1070x10-3 
-0.3709 -0.0205 0.5412 x10-3 
3.1 脉冲 激励 响应 
假定 在 搁 性 梁 的 自由 端 端点 X 方 向 受到 脉冲 力 


和 矩 激励 ,仿真 计算 该 端点 工 方 向 和 了 方向 的 位 移 输 
上 响应 。 取 观测 矩阵 Cs 和 C, .控制 力 分 布 矩阵 瑟 、 
输出 了 中 各 分 量 相 对 重要 性 的 权 和 矩阵 2 和 阻尼 比 < 
分 别 为 
Ci=03x16, C= [03 .3 : 
=[0 0 0 : 1 
QO=1st =0,.01 
根据 数值 解 式 (16) 计 算得 到 模 态 价值 ,根据 式 
(19) 计 算得 到 MQI 值 ,结果 如 表 4 所 示 , 并 给 出 其 
柔性 振动 的 前 10 阶 模 态 作为 梁 原 系统 的 真实 响应 。 


于 


1 :Qi0d 和 
0 0 0 .oo] ， 
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表 4 脉冲 激励 下 模 态 价值 计算 结果 


Tab.4 Results of modal cost analysis under impulse excitation 


学 报 


i wi/Hz w/ Hz V MOQI 

1 0.19345 0.2085 2.6481x10-4 0.9973 

2 1.3781 1.3906 5.7619x1077 0.0022 

3 1.6479 1.6494 1.1629x10- 4.3794x10-8 
4 4.0587 4.0878 1.3504x10-’ 5.0855x10-4 
5 5.4160 5.4180 1.4914x10-! 5.6165x10-° 
6 8. 162 4 8.1709 1.4503x10-” 5.4618x10 
7 9.9986 9.9993 ”4.5741x10-2 1.7226x10-8 
8 13.7077 13.7621 1.3648x10-® 5.1398x10-5 
9 15.4739 15.4770 2.5793x10- 9.7135x10-8 
10 20.7316 20.7205 4.2334x10-' 1.5943x10- 


由 表 4 可 知 ,系统 最 大 的 模 态 价值 是 第 1 阶 , 其 
MQI =0.997 3 ,达到 总 价值 的 99.73% ,对 系统 贡献 最 
大 ,其 他 阶 模 态 与 第 1 阶 模 态 相 比 MQI 数值 上 有 数量 
级 的 差异 。 故 选择 第 1 阶 模 态 为 模 态 价值 降 阶 系统 。 

同时 ,利用 方程 式 (15 ) 进行 内 平衡 降 阶 ,计算 
导 内 平衡 系统 的 奇异 值 为 
o=[0.1782 0.1782 0.0032 0.0032 
0.0009 0.0009 0.0002 0.0002 
6.6442x10” 6.6433x10” 2.2557x10™ 
1.3125x10™” 1.3063x10™”” 6.2664x10™” 6.2447x10™ 
3.3285x10 3.3263x10 2.5973x10™* 2.5946x10™] 

由 上 述 计算 得 到 的 奇异 值 可 知 ,第 1 阶 模 态 对 
应 的 两 个 奇异 值 最 大 ,其 他 阶 模 态 与 第 1 阶 模 态 相 
比 数值 上 有 数量 级 的 差异 。 故 取 第 1 阶 内 平衡 模 态 
为 内 平衡 降 阶 系统 。 

取 输 入 为 脉冲 信号 ,脉冲 强度 为 10, 周 期 为 
0. 01。 利 用 非 约 束 模 态 平 动 耦合 系数 、 转 动 耦合 系 
数 和 降 阶 准则 进行 数值 仿真 。 原 系统 和 降 阶 系统 的 
响应 如 图 2 ~4 所 示 。 从 结果 看 出 , 降 阶 系统 与 原 系 
统 的 响应 吻合 良好 。 


0.015 


2.2473 x10™ 


i Wp 


本 
po et 


i 昌 
0 
nl 


| 和 nu 
ro 
时 


响应 /m 
> 
Ey 
3 


一 一 方向 
| 


500 1 000 
时 间 /s 
图 2 脉冲 激励 下 原 系统 挠 性 梁 未 端 响应 


Fig.2 Tip response of the flexible beam by using the original 


1500 2 000 


system under impulse excitation 
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汪 
0.015 表 5 正弦 激励 下 模 态 价值 计算 结果 


Tab.5 Results of modal cost analysis 


under sinusoidal excitation 


Bd 
Cd 


wi/Hz ws/ Hz V MQI 


上 0 tt 1 0.1935 0.2085 1.3127x10-3 1.1330x1075 
加 0.000 pm uh 
可 hs 2 1.3781 1.3906 3.5861x10-! 3.0952x10- 
-0.005 ， 3 1.6479 1.6494 1.0997x10-8 0.949 2 
shh a 4 4.0587 4.0878 3.6580x10-4 3.1572x10-6 
R000 5 5.4160 5.4180 5.8392x10-0 0.0504 
6 8.1624 8.1709 2.8256xl10-4 2.4388x10-5 
ee qh S090 7 9.9986 9.9993 1.3184x10-3 1.1379x10-5 
Ns 8 13.7077 13.7621 2.9514x10-16 2.5474x10- 
图 3 脉冲 激励 下 模 态 价值 降 阶 系统 挠 性 梁 未 端 响应 9 15.4739 15.4770 1.2790x10-? 1.1039x1074 
Fig.3 Tip response of the flexible beam by using the modal 10 20.7316 20.7205 13914x10- 1.2009 x107° 
cost analysis reduction system under impulse excitation 由 上 述 计 算得 到 的 奇异 值 可 知 ,第 1 阶 模 态 对 
0.015 应 的 两 个 奇异 值 最 大 , 故 取 第 1 阶 内 平衡 模 态 为 内 
™ ee > 
my 平衡 降 阶 系统 。 
GN ss 取 输 入 为 正弦 外 激励 , 其 振幅 为 2, 频率 为 
0.005 区 向 和 0.010 ， 4 
0 i i Bh 1 rad/s。 仿 真 所 得 到 的 原 系 统 响 应 和 降 阶 系统 响应 
总 hi i! Wt MFT = 
Co 好 如 图 5 -7 所 示 。 根 据 图 5 ~ 7 可 知 , 降 阶 系统 与 原 
全 机 TH > Se EN 
= 和 系统 吻合 良好 , 故 两 种 准则 均 能 够 起 到 较 好 的 降 阶 
NN 一 坊 册 效果 。 
= -0.010 : 
CN 
-015 
NN 500 1 000 1 500 2 000 
2 时 间 /s 


二 图 4 脉冲 激励 下 内 平衡 降 阶 系统 挠 性 梁 末 端 响应 
Eiga=4 Tip response of the flexible beam by using the internal 


响应 /(10 my) 


nce technique reduction system under impulse excitation 


3, 守 正弦 激励 响应 


Ra 2 
@ 假 设 在 挠 性 梁 的 自由 端 端 点 了 方向 受到 正弦 力 NA 
矩 激励 ,仿真 计算 该 端点 蕊 方向 和 了 方向 的 位 移 答 到 5 正弦 激励 下 原 系统 挠 性 染 未 端 响应 
出 响应 。 取 控 制 力 分 布 矩 阵 B 如 下 ,其 他 数量 保持 Fig.5 Tip response of the flexible beam by using the 
不 变 , 即 original system under sinusoidal excitation 
~ 下 ， 


B=[0 0 0 : 0 10 : 02o] 
表 5 采用 数值 解 式 (16) 计 算得 模 态 价值 结果 
从 表 5 中 可 以 看 出 , 模 态 价值 最 大 值 出 现在 第 3 
阶 模 态 ,其 MQI=0.9492, 对 系统 贡献 最 大 。 故 选择 
第 3 阶 模 态 为 模 态 价值 降 阶 系统 。 
计算 得 内 平衡 系统 奇异 值 为 
or=[0.0004 0.0004 5$.188$x10 5.1853x10™ 
7.9352 x10- 7.9016x10-® 3.9672x10- 3.9668 x10-° 


响应 /(10 m) 


时 间 /s 
1.43$8x10 1.4336x10™® 5.7086x10-” 5.6974x107 图 6， 正 蓄 激 励 下 模 态 价值 降 阶 系统 挠 性 梁 未 端 响应 
4.7038x10-” 4.6915x10-” 2.9154x10-” 2.9074x10™ Fig.6 Tip response of the flexible beam by using the modal 
1.2916x10-” 1.2845x10-” 2.2369x10™ 2.2353x107] cost analysis reduction system under sinusoidal excitation 
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响应 /(10 mm) 


时 间 /s 
图 7 正弦 激励 下 内 平衡 降 阶 系统 挠 性 梁 末 端 响应 
Fig.7 Tip response of the flexible beam by using the internal 


balance technique reduction system under sinusoidal excitation 
疆 -从 
4 纺 论 


全 本 研究 针对 中 心 刚体 - 挠 性 梁 模 型 的 挠 性 航天 
过 非 约 束 模 态 振 型 展开 ,得 到 非 约 束 模 态 描述 
系统 刚体 运动 和 挠 性 变形 的 动力 学 方程 。 结 合 模 态 
侍 分 析 准则 和 内 平衡 降 阶 准则 分 析 了 挠 性 航天 器 
模型 降 阶 问题 。 模 态 价值 分 析 准则 讨论 了 各 阶 模 态 
的 模 态 价值 对 系统 贡献 的 重要 程度 , 因此 可 以 保留 
Fb 染 统 贡献 较 多 的 模 态 ,从 而 组 成 降 阶 系统 。 内 平 
f 降 阶 准则 通过 非 奇 异 线 性 变换 ,将 状态 空间 系统 
畦 他 为 内 平衡 系统 ,从 而 保留 该 系统 内 可 控 可 观 
Gram 对 角 矩 阵 对 角 线 上 较 大 的 值 所 对 应 的 模 态 来 
对 又 统 进行 模 态 截断 ,从 而 达到 降 阶 效果 。 最 后 通 
过 和 蒋 真 实验 验证 了 两 个 准则 对 非 约束 模 态 动力 学 模 
本 用 认 与 和 的 有 效 性 ,结果 表明 模 态 从 人 分 析 兴 
i 内 平衡 降 阶 准则 也 能 对 非 约束 模 态 动力 学 模型 
ns 
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